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Abstract: The third del Pezzo surface admits a unique Kahler-Einstein metric, which is not 
known in closed form. The manifold's toric structure reduces the Einstein equation to a single 
Monge- Ampere equation in two real dimensions. We numerically solve this nonlinear PDE 
using three different algorithms, and describe the resulting metric. The first two algorithms 
involve simulation of Ricci flow, in complex and symplectic coordinates respectively. The 
third algorithm involves turning the PDE into an optimization problem on a certain space 
of metrics, which are symplectic analogues of the "algebraic" metrics used in numerical work 
on Calabi-Yau manifolds. Our algorithms should be applicable to general toric manifolds. 
Using our metric, we compute various geometric quantities of interest, including Laplacian 
eigenvalues and a harmonic (1, l)-form. The metric and (1, l)-form can be used to construct 
a Klebanov-Tseytlin-like supergravity solution. 
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1. Introduction 

Kahler metrics on manifolds play an important role in mathematics and physics. As Yau 
demonstrated [1], in the Kahler case it is often possible to prove the existence (or nonexistence) 
of metrics which solve the Einstein equation. While it is extremely valuable to know whether 
they exist, for many purposes one also wants to know their specific form. The existence 
theorems, however, are generally non-constructive, and explicit examples of Kahler-Einstein 
metrics are rare. This state of affairs naturally leads to the following question: Is it possible, 
in practice, to find accurate numerical approximations to these metrics using computers? The 
last two years have seen significant success, with a variety of different algorithms providing 
numerical solutions to the Einstein equation on the Calabi-Yau surface K3 [2-4] and on a 
three-fold [5]. The Kahler property proved to be as crucial for this numerical work as it was 
for the existence theorems. 

In this paper we extend this success to a non-Calabi-Yau manifold, namely CP^ blown up 
at three points (CP^#3CP ), also known as the third del Pezzo surface (dP^). This manifold 
is known by work of Siu [7] and Tian-Yau [6] to admit a Kahler-Einstein metric with positive 
cosmological constant, but as in the Calabi-Yau case that metric is not known explicitly. An 
important property of dP^ that differentiates it from Calabi-Yau manifolds is that it is toric. 
Toric manifolds are a special class of Kahler manifolds whose ?7(1)"' isometry group (where n 
is the manifold's complex dimension) allows even greater analytical control, and we develop 
algorithms for solving the Einstein equation that specifically exploit this structure. We should 
also note that the Kahler-Einstein metric on dP^ is unique (up to rescaling); therefore, rather 
than having a moduli space of metrics as we have in the case of Calabi-Yau manifolds, there 
is only one metric to compute. 

On the physics side, the Kahler-Einstein metric on dP^ is important because it can be 
used to construct an example of a gauge/gravity duality. These dualities provide a bridge 
between physical theories of radically different character, allowing computation in one theory 
using the methods of its dual. Given the Kahler-Einstein metric on dP^ , one can construct a 
five-dimensional Sasaki-Einstein metric. Compactifying type IIB supergravity on this man- 
ifold one obtains an AdSs supergravity solution, which has a known superconformal gauge 
theory dual [8-10]. An interesting generalization includes a 3- form flux from wrapped D5- 
branes. This flux can be written in terms of a harmonic (1, l)-form on dP^, which we also 
compute numerically in this paper. The resulting supergravity solution (the analogue of the 
Klebanov-Tseytlin solution on the conifold [11]) is nakedly singular, but is the flrst step to- 
ward flnding the full supergravity solution on the smoothed-out cone (the analogue of the 
Klebanov-Strassler warped deformed conifold [12]). This smoothed-out cone is dual to a 
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cascading gauge theory, and knowing the exphcit form of the supergravity solution would 
be useful for both gauge theory and cosmology applications. This physics background is 
explained in detail in Section |2|. 

In Section ^ we review the mathematical background necessary for understanding the 
rest of the paper. Here we closely follow the review article on toric geometry by Abreu [13]. 
We explain the two natural coordinate systems on a toric manifold, namely complex and 
symplectic coordinates. In complex coordinates, the metric is encoded in the Kahler potential, 
and in symplectic coordinates in the symplectic potential; these two functions are related 
by a Legendre transform. In either coordinate system, the Einstein equation reduces to a 
single nonlinear partial differential equation, of Monge-Ampere type, for the corresponding 
potential. Thanks to the U{1)"' symmetry, this PDE is in half the number of dimensions of 
the original manifold (two real dimensions for dP^). In this section we also derive the equation 
we need for the (1, l)-form, and give all the necessary details about dP^. 

In this paper we describe three different methods to solve the Monge-Ampere equation. 
In Section § we explain the first two methods, which involve numerically simulating Ricci flow 
in complex and symplectic coordinates respectively. Specifically, we use a variant of Ricci flow 
(normalized) that includes a A term, 

^ = -2R^, + 2Ag^,, (1.1) 

whose fixed points are clearly Einstein metrics with cosmological constant A. According to 
a recent result of Tian-Zhu [14], on a manifold which admits a Kahler-Einstein metric, the 
flow (^]^) converges to it starting from any metric in the same Kahler class. Our simulations 
behaved accordingly, yielding Kahler-Einstein metrics accurate (at the highest resolutions 
we employed) to a few parts in 10^ (and which agree with each other to within that error). 
Numerical simulations of Ricci flow have been studied before in a variety of contexts [15-17], 
but as far as we know this is the first time they have been used to find a new solution to the 
Einstein equation (aside from a limited exploration of its use on K3 [2]). In this section we 
also explore the geometry of this solution, and discuss a method for computing eigenvalues 
and eigenfunctions of the Laplacian in that background. 

In Section ^ we introduce a different way to represent the metric, based on certain poly- 
nomials in the symplectic coordinates. This non-local representation is a symplectic analogue 
of the "algebraic" metrics on Calabi-Yau manifolds employed in numerical work by Donald- 
son [3] and Douglas et al. [4,5]. To demonstrate the utility of this representation we give a 
low order polynomial fit to the numerical solutions found by Ricci flow, that can be written 
on one line and yet agrees with the true solution to one part in 10^, and everywhere satisfies 
the Einstein condition to better than 10%. 

In Section ^ we discuss our third method to compute the Kahler-Einstein metric. As 
above we represent the metric using polynomials in the symplectic coordinates, but now we 
constrain the polynomial coefficients by solving the Monge-Ampere equation order by order in 
the coordinates. This leaves a small number of undetermined coefficients, which we compute 
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by minimizing an error function. Using this method we obtain numerical metrics of similar 
accuracy to those found by Ricci flow. By a similar method we also calculate eigenvalues and 
eigenfunctions of the Laplacian, and the harmonic (1, l)-form. 

We conclude the paper with a brief discussion of the three methods and their relative 
merits in Section |^. 

At the two websites [21], we have made available for download the full numerical data 
representing our metrics, as well as Mathematica notebooks that input the data and allow 
the user to work with those metrics. The codes used to generate the data are also available 
on those websites. 

We believe that all of the methods we present here can be applied to a general toric 
manifold. We will report elsewhere on an application to dP2 [18], which does not admit a 
Kahler-Einstein metric but does admit a Kahler- Ricci soliton [19]. For future work, there 
is also a natural analogue of dPa to study in three complex dimensions. From Batyrev's 
classification of toric Fano threefolds, it follows that there are precisely two of them that admit 
Kahler-Einstein metrics which are not themselves products of lower dimensional manifolds [20, 
Section 4]. One of these is CP^. The other is the total space of the projectivization of the 
rank two bundle O 0(1, -1) over CP^ x CP^ Its Delzant polytope (see Section | for the 
definition) possesses a D4, symmetry, and a fundamental region is simply a tetrahedron. 

While this work was in progress we have learned that Kahler-Einstein metrics on dP^ 
have also been computed in [22] using the methods of [3]. 

2. Gauge/gravity duality 

A prototypical example of a gauge/gravity duality which provides the physics motivation 
for studying dP-^ is the Klebanov-Strassler (KS) supergravity solution [12]. (Mathematicians 
may wish to skip this section.) The KS solution is a solution of the type IIB supergravity 
equations of motion. The space-time is a warped product of Minkowski space M^'^ and the 
deformed conifold X. The affine variety X has an embedding in defined by 

jZ^l = e (2.1) 

i=l 

where S C. There are also a variety of nontrivial fluxes in this solution which we will return 
to later. 

One important aspect of the KS solution is its conjectured duality to a non-abelian gauge 
theory, namely the cascading SU{N) x SU{N + M) M = 1 supersymmetric gauge theory 
with bifundamental fields Ai and Bi, i = 1 or 2, and superpotential 

W = e,jekiA,BkAjBi . (2.2) 

This theory is similar in a number of respects to QCD; it exhibits renormalization group 
flow, chiral symmetry breaking, and conflnement. Moreover, all of these properties can be 
understood from the dual gravitational perspective. 
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In addition to its gauge theory applications, the KS solution is important for cosmology. 
Treating the deformed conifold as a local feature of a compact Calabi-Yau manifold, the KS 
solution provides a string compactification with a natural hierarchy of scales in which all the 
complex moduli are fixed [23]. Stabilizing the Kahler moduli as well [24], the KS solution 
can become a metastable string vacuum and thus a model of the real world. In this context, 
inflation might correspond to the motion of D-branes [25] and cosmic strings might be the 
fundamental and D-strings of type IIB string theory [26]. 

One naturally wonders to what extent the cosmological and gauge theoretic applications 
depend upon the choice of the deformed conifold. A natural way to generalize X is to consider 
smoothings of other Calabi-Yau singularities. One such family of singularities involves a 
Calabi-Yau where a del Pczzo surface dPn shrinks to zero size. (Here dPn is CP^ blown up 
at n points.) Note we are distinguishing here between resolutions — or Kahler structure 
deformations — where even dimensional cycles are made to be of finite size, and smoothings 
— complex structure deformations — where a three dimensional cycle is made finite. Using 
toric geometry techniques, Altmann [27] has shown that the total space of the canonical 
bundle over dPi admits no smoothings, while dP2 admits one and dP'^ two.^ The higher 
dPn are not toric. Thus two relatively simple candidates for generalizing the KS solution are 
smoothed cones over the complex surfaces dP2 and dP^. 

Without knowing the details of the metric on the smoothed cone X, one can show that 
a generalization of the KS solution exists for such warped products [32]. The solution will 
have a ten dimensional line element of the form 

ds^ = h{py^''^r]f,ydx^'dx'' + hipf^dslc , (2.3) 

where ds\ is the line element on X, p G X, the map h : X ^ is called the "warp factor", 
and rjni, is the Minkowski tensor for M}'^ with signature ( — h ++). There are also a variety 
of nontrivial fluxes turned on in this solution. There is a five form flux 

F5 = dCi + ★dC4 , where C4 = -^dx^ A dx^ A dx'^ A dx^ , (2.4) 

and where gs is the string coupling constant. The finite smoothing indicates the presence of 
a harmonic (2, l)-form lj2,i which we take to be imaginary self-dual: *x'^2,i = ^"^2,1- From 
^2,1 we construct a three-form flux G3 = Cuj2,i where C is a constant related to the rank of 
the gauge group in the dual theory. The warp factor satisfies the relation 

l^xh = -^^GaUG*f\ (2.5) 

where the indices on G3 are raised and the Laplacian Ax is constructed using the line element 
ds\ without the warp factor. Although this solution holds for general X, clearly in order to 

^Tho physical relevance of this fact for supersymmetry breaking was pointed out in [28-30] (see also [31] 
for a more recent account). 
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know detailed behavior of the fluxes and warp factor as a function of p, we need to know a 
metric on X. 

Prom this perspective, we have chosen to study dP3 and not dP2 in this paper because 
dP^ is known to have a Kahlcr-Einstein metric [6] while dP2 does not [33]. Given that dP^ is 
Kahler-Einstein, we can construct a singular Calabi-Yau cone over dPs in a straightforward 
manner: 

ds\ = dr^ + [{di; + af + dsl] , (2.6) 

where a = —2i{df — df) and / is half the Kahler potential on V = dP^. Here dsy is a 
Kahler-Einstein line element on dPs. In a hopefully obvious notation, r is the radius of the 
cone and an angle. Although such a cone over dP2 probably exists as well, it will involve an 
irregular Sasaki-Einstein manifold as an intermediate step; the metric on the Sasaki-Einstein 
manifold over dP2 is not yet known. 

Given that dP^ is Kahlcr-Einstein, the problem of finding aj2,i on the singular cone 
reduces to finding a harmonic (1, l)-form 6 on dP^ such that 6 A lv = (where ui is the 
Kahler form on V) and -kyS = —Q, as pointed out in [34]. The relation between a;2,i and 9 is 
0^2,1 = {—idr/r + dtp + a) AO. 

In addition to finding a numerical Kahler-Einstein metric on dP^, we will also find a 
numerical (1, l)-form 9, thus yielding a singular generalization of the KS solution for dP^. 
Historically, before the KS solution, Klebanov and Tseytlin [11] derived exactly such a singular 
solution for the singular conifold. Although we have not produced a numerical solution for 
h{p), with the explicit metric and numerical (l,l)-form for dP^ in hand, we have all the 
necessary ingredients to calculate h{p). The KT solution for dPi is known [35]. 

The next step would be to find a Ricci flat metric and imaginary self-dual (2,l)-form 
on the smoothed cone over dP^, thus providing a generalization of the KS solution. Such a 
solution would open up many future directions of study, both in gauge theory and cosmology. 
To name a handful of possibilities, one could compute A;-string tensions of the confining low 
energy gauge theory dual to this dPs background, generalizing work of [36]. Alternatively, 
treating the SUGRA solution as a cosmology, one could compute annihilation cross sections 
of cosmic strings [26] or slow roll parameters for D-brane inflation [25]. 

3. Mathematical background 

3.1 Complex and symplectic coordinates 

The formalism we use to construct our numerical Kahler-Einstein metric on dP^ is based on 
work by Guillemin [37] and later developed by Abreu [13]. Here we summarize this formalism. 

We consider an n (complex) dimensional compact Kahler manifold M. The manifold is 
equipped with the following three tensors: 

• A complex structure J^^i, satisfying J^^^J'^x = ^<^^A- 

• A symplectic form (also called in this context a Kahler form) lo, which is a non- 
degenerate closed two-form. 
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• A positive-definite metric ds^ = g^^dx^dx'^ . 

These tensors are related to each other by: 

3iiv = ^ixxJ^v (3.1) 

Now let M also be a toric manifold. Tensors which are invariant under its U (1)" = T" 
group of diffeomorphisms we call toric. In particular, we will restrict our attention to toric 
metrics. Let M° be the subset of M which is acted on freely by that group. There is a natural 
set of n complex coordinates z = u + iO on M°, where u G M" and 9 £ T"^ {6i 9i + 2it); the 
?7(1)" acts on 9 and leaves u fixed. Given that the manifold is Kahler, the metric may locally 
be expressed in terms of the Hessian of a Kahler potential f{z)? Because g^^i, is invariant 
under the action of C/(l)", the potential can be chosen to be a function of u. The line element 
is 

ds^ = gijdzidzj + gydzidzj = Fij [duiduj + d9id9j), (3-2) 
where we have introduced Fij{u): 

_ _ 1 _ 1 

The Kahler form is 

u = 2iddf = igijdzi A dzj = dui A dOj . (3-4) 
The complex structure in these coordinates is trivial: 

= J%=5i, J\ = J% = 0. (3.5) 

It is often convenient to work with symplectic coordinates w = x + iO, which are related 
to the complex coordinates by: 

X = - . (3.6) 

Under this map (also known as the moment map), is mapped to the interior P° of a 
convex polytope P C which is given by the intersection of a set of linear inequalities, 

P = {x: la{x) > Va} , la{x) =Va-X + Xa] (3.7) 

the index a labels the faces, and the normal vector Va to each face is a primitive element 
of Z". These Va define the toric fan in the complex coordinates z. We will see that the Aa 
determine the Kahler class of the metric. In terms of these new coordinates, the Kahler form 
becomes trivial: 

uj = dxiA dOi . (3.8) 
Introducing the symplectic potential, which is the Legendre transform of /, 

g{x) = u-x- f{u) , (3.9) 



^We follow the conventions of Abreu [13]. Note that / is one-half the usual definition of the Kahler potential. 
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the line element can be written 



ds^ = Gij dxidxj + G'^ dOiddj , (3.10) 

where 



Cidx-j 

and G*-' is the inverse of Gij. Note that, while F^j is regarded as a function of u and G*-' 
as a function of x, under the mapping ( |3.6| ) the two matrices are equal to each other. The 
complex structure in symplectic coordinates is given by: 

r^g^=-G'^, j'\^=Gij, J\ = J\ = 0. (3.12) 

To summarize, the complex and symplectic coordinate systems are related to each other 
by a Legendre transform: 



du ' dx 



f{u)+g{x)=u-x, F,,{u) = G'^{x). (3.13) 



Note that the whole system (complex and symplectic coordinate systems) has four gauge 
invariances under which the metric is invariant and the relations ( |3.13| ) are preserved: 

1. /(n) — > f{u) + c, g{x) g{x) — c for any constant c; 

2. /(n) — > f{u) + /c • u, X — > x + A; for any vector k; 

3. g{x) — > g{x) + k ■ x, u ^ u + k ioi any vector k; 

4. x — > M • x, li* — > li* • M"-*^ for any element M G GL(n, Z) (the ring Z is necessary to 
preserve the integrality of the boundary vectors Va)- 

However, gauge invariances (2) and (4) are broken by the polytope down to the subgroup of 
M" X GL(n, Z) under which P is invariant. 

3.2 Boundary conditions and the canonical metric 

The complement of M° in M consists of points where one or more circles in the fiber 
degenerate. In order to have a smooth metric on all of M, there are two boundary conditions 
that must be imposed on f{u), or equivalently g{x). These are somewhat easier to express 
in the symplectic coordinate system. The first condition is that, as we approach a face of 
the polytope, the part of the metric parallel to the face, should not degenerate (or become 
infinite). Technically the requirement is that the function 



det Gij Y[la, (3.14) 



which is positive in the interior of the polytope, should extend to a smooth positive function 
on the entire polytope. The second condition is that the shrinking circle(s) should go to zero 
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size at the correct rate in order to avoid having a conical singularity. To express this boundary 
condition in terms of g, Guillemin [37] and Abreu [13] introduced the canonical symplectic 
potential: 

9can= -^^la^T^la- (3.15) 

a 

This canonical potential leads to a metric on M that is free of conical singularities. Further- 
more, every smooth metric corresponds to a 5 that differs from g^an by a function that is 
smooth on the entire polytope; we will call this function h: 

9 = 5can + h. (3.16) 

So far, we considered a toric manifold with a fixed metric. In general, a given toric 
manifold will admit many different toric metrics, i.e. Kahler metrics invariant under the 
given f/(l)" diffcomorphism group. These will be described by functions g(x) (or f{u)) that 
differ by more than the gauge transformations listed above. Metrics with the same polytope 
(modulo the gauge transformations 2 and 4, which act on the polytope) are in the same Kahler 
class. In general, the topology of the manifold is determined by the number and angles of the 
polytope's faces, i.e. the vectors Va (modulo gauge transformation 4), while the Kahler class 
is determined by their positions, i.e. the numbers Xa (modulo gauge transformation 2). If 
the Xa are all equal, then the Kahler class is proportional to the manifold's first Chern class. 
More specifically, if Xa = for all a, then A[a;] = 27rci(M). The case of interest, dPs, has 
ci(M) > 0, so we must take A > 0. 

3.3 Examples 

In one complex dimension, there is only one compact toric manifold, CP^. The corresponding 
polytope P is simply the interval, whose length determines the Kahler modulus. We take 
P = [—A, A]. The canonical symplectic potential 

5can(x) = ^ ((A + x) ln(A + x) + {X-x) ln(A - x)) (3.17) 

yields the round metric of radius \/A: 

A^ - 

dScan = T2 o^^"^ + ' (^"l^) 

A — X A 
The Kahler coordinate u is related to x by 

X = Atanhu, (3.19) 



and the Kahler potential is 
giving the metric in the form 



fca.n{u) = A In cosh u, (3.20) 
c^Scan = Asech^ u{du^ + de'^). (3.21) 
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Figure 1: The polytopes for the five compact toric manifolds with positive first Chern class: from 
left to right, CP^ x CP\ CP^ dPi, dP2, dP^. Each polytope is drawn such that the for all its faces 
are equal, corresponding to the Kahler class being proportional to the first Chern class. 



There are five compact toric surfaces with positive first Chern class (i.e. toric Fano sur- 
faces); their polytopes are shown in Fig. ||. CP""^ x CP^ has two moduh, the sizes of the two CP^ 
factors; when these are equal the canonical metric is Einstein. CP^ has only a size modulus, 
and again the canonical metric is Einstein, as shown in Appendix The del Pezzo surfaces 
dPi, dP2, and dP^j have two, three, and four moduli respectively. Their canonical metrics 
are never Einstein. Indeed, dPi and dP2 do not admit Kahler-Einstein metrics at all [33], 
essentially because their Lie algebra of holomorphic one-forms is not reductive. The case of 
dP^ is special. On the one hand, it is known to admit a toric Kahler-Einstein metric [6]. On 
the other hand, that metric is not the canonical one; indeed it is not known in closed form 
— hence the necessity of computing it numerically. dP^ is discussed in more detail in Section 
below. 



3.4 The Monge- Ampere equation 

As usual for a Kahler manifold, the Ricci curvature tensor Rij can be written in complex 
coordinates z in a simple way, namely 

Recall that there is a two-form TZ = iRijdz'' A dP associated with Rij where the class [TZ] = 
2ttci{M). We are interested in Kahler-Einstein metrics on M, that is metrics which satisfy 
the following relation 

Rij = Agij (3.23) 

for some fixed A. As explained above, this implies that = A~^ for all a. The sign of A is 
not arbitrary but is fixed by the first Chern class of M, which we are now assuming to be 
positive. 



Given ( |3.22|) , we may integrate ( |3.23|) twice, yielding 

In det Fij = -2A/ + 7 • -u - c (3.24) 
where 7 and c are integration constants. In symplectic coordinates, this becomes 

In det Gii = -2Ag + ^ • (2Ax - 7) + c . (3.25) 

ox 
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Both ( p. 241) and ( |3.25| ) are examples of Monge- Ampere type equations, which must be solved 
with the boundary conditions discussed in Section ^]2| above. The values of the constants 7 
and c in these equations are arbitrary; given a solution with one set of values, a solution with 
any other set can be obtained using gauge transformations (1) and (2). However, as discussed 
above, gauge transformation (2) acts on the polytope P by a translation. Therefore, if we fix 
the position of P, there will be a unique 7 such that ( 3.25| ) admits a solution. 



Since the boundary conditions are non-standard, it is worth exploring them in more 
detail. Again, we work in the symplectic coordinate system. Recall that the condition for a 
smooth metric was that h{x), defined by 

g = S-can + h, (3.26) 
be smooth on P, including on its boundary. We can re-write ( |3.25| ) in terms of h as follows: 

Indet + Gt^^^ ] = -2Ah + - • (2Ax - 7) - Pcan + c , (3.27) 



Pcan = In det G^f + 2A5can " ^ " (2Ax - 7) . (3.28) 



where 

dgc 
dx 

Normally, for a second-order PDE, we would expect to have to impose, for example, Dirichlet 
or Neumann boundary conditions in order to obtain a unique solution. However, in the case of 
( |3.25| ), the coefficient of the normal second derivative goes to zero (linearly) on the boundary 
(near the face a of the polytope, G^^^v^ ~ la)- Therefore, under the assumption that h and 
its derivatives remain finite on the boundary, the equation itself imposes a certain (mixed 
Dirichlet-Neumann) boundary condition. If we were to try to impose an extra one, we would 
fail to find a solution. This is illustrated by the case of CP^ Setting A = A-^ we have 



Pcan = 0, so that ( 3.27 ) becomes 



In (^1 + ^^-J^h"^ = ji-h + xh') + c; (3.29) 

we see that the coefficient of h" vanishes on the polytope boundary. 

For future reference we record here the formulas for the Ricci and Riemann tensors in 
symplectic coordinates. Their non-zero components are 



^XiXjX]^xi — ^XiXj — -ft xj^xi — -fi- — '^ml[i'^j]nk {o.ol) 



TjOi Ok — _ I? ^j"k — _ Tf^i — P t'j t'l 
J L -y. - -y. 7 J l/T - -y. 7 J L ■ T" 7 (^T* - T* 7 

— 2 {'^ijkl ~ ^ '^ijm^kln — ^ ^ml{i^j)nk) ; {o.OZ) 

where we've defined 
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3.5 Volumes 

Here follows a short discussion about volumes useful for understanding the relation between 
Xa and A above. 

We know that the volume of a complex surface M is 

Vol(M) = Lo\ (3.34) 

while the volume of a curve C is 

Vol(C7) = [ io . (3.35) 
Jc 

From the Kahler-Einstein condition ( p. 23] ), ( p. 41 ), and the fact that the class of the Ricci form 
is related to the first Chern class, [TZ] = 2ttci{M), it follows that [uj] = 27rci(M)/A and that 

Vol(M) = =^ci{Mf . (3.36) 

For CP^ blown up at k points, cf = 9 — k. Meanwhile, for our curve, 

27r 

Vol(C) = — ci(M) • C7 . (3.37) 



In symplectic coordinates, it is easy to compute these volumes. From (3^) or ( p.lOD , 
the volume (p.34|) reduces to 47r^ times the area of P. Setting A = 1 corresponds to setting 



the Xa = 1. For curves, the computation is similarly easy. For example, some simple torus 
invariant curves correspond to edges of P, and the volume computation reduces to measuring 
the length of an edge of P. 

3.6 The Laplacian 

The Laplacian acting on a scalar function ip is, 

1 



AtA = rr— V^^^tg^ g'^'d.i; . (3.38) 



In symplectic coordinates, the determinant of the metric y^detg^ = 1. We consider a 
function ■0 invariant under the torus action, implying no dependence on the two angular 
coordinates so that ^/J = 'tp{xi,X2)- Thus, the Laplacian can be written in a simpler fashion: 

The Laplacian thus depends on a third derivative of g. However since we are interested in 



AeV' = G'i^^ - 2Ax,^ . (3.40) 



Einstein metrics, we may take ( |3.25| ), differentiate it, and use it to eliminate this derivative. 
One then finds the form 

_^!^ _2Ax- — 
dxidxj ^ dxi 

From this equation one deduces the interesting fact that the symplectic coordinates are eigen- 
functions of — Ae with eigenvalue 2 A — the symplectic coordinates are close to being harmonic 
coordinates. 
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3.7 Harmonic (1, l)-forms 

A compact toric variety M of dimension n defined by a fan {va} of m rays has Betti number 
1)2 = m — n [38]. Moreover, dPa has no holomorphic (nor antiholomorphic) 2-forms and thus 
must have m — n harmonic (l,l)-forms. These (l,l)-forms are in one-to-one correspondence 
with torus invariant Weil divisors Da modulo linear equivalence. The equivalence relations 
are 

(3.41) 

a 

From their connection with the divisors Da, perhaps it is not surprising that these (1,1)- 
forms 0a can be constructed from functions which have a singularity of the form ln{va-x+l) 
along the boundaries of the polytope [13]. Moreover, they should satisfy Maxwell's equations, 

^6* = and d*6' = . (3.42) 

The first equation d9 = is immediate from the local description of Oij as didjfi. The other 
equation is 

= D'Oij = g'^'D-^Oij . (3.43) 
Using the Bianchi identity, we find then 

= g''Dj9^-, = dj [g%j,) , (3.44) 

or 

g^'^Oij. = constant . (3.45) 

Now this last equation is rather interesting. The left hand side is the Laplacian operator 
acting on /i. In symplectic coordinates, the left hand side can be rewritten to yield 

^ (g'^-^] = constant . (3.46) 

OXi \ OXj J 



Using the Kahler-Einstein condition, the left hand side becomes ( 3.4C1| ). 
The constant is easy to establish. We have that 

Vol(M)5*%-= /" 9A*uj. (3.47) 

The Kahler form is self-dual under the Hodge star, ★a; = to. Conventionally, we may write 
[9] = 2ttci{D), assuming 9 is the curvature of a line bundle 0{D). Using finally that [uj] = 
27rci(M) (for A = 1), we find that 

^?% = 2^, (3.48) 
where K is the canonical class of M. This constant is often referred to as the slope of 0{D). 
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3.8 More on dP^ 

For dPs, the toric fan is described by the six rays spanned by 

VI = {1,0); 1-2 = (1,1); vs = {0,l); 



(3.49) 



V4 



(-1,0); V5 = i-1,-1); ve = iO,-l). 



This fan leads to a dual polytope P which is a hexagon. As discussed above, when all six 
Xa are equal, the Kahler class is proportional to the first Chern class; we will choose Aa = 1. 
In this case the polytope has a dihedral symmetry group Dq, generated by the following Z2 
reflection and rotation: 



1 

1 



R2 



1 1 

-1 



(3.50) 



This discrete symmetry will be shared by the Kahler-Einstein metric on dP^, and will therefore 
play an important role in our computations. Note that the element i?2 acts which 
sets 7 = in ( ^.241 ) and ( 3.25 ). Note also that the acts naturally on the group of Cartier 
divisors on dP^ (and hence on the harmonic (l,l)-forms), as can easily be seen by thinking 
of the set of Va as divisors on the manifold. 

The hexagon P has a natural interpretation as the intersection of the polytope for a 
symmetric unit (CP"^)^, which is a cube with coordinates xi,X2,X3 satisfying \xi\ < 1, with 
the plane xi + X2 + x^ = 0. The intersection defines a natural embedding of dP^ into (CP^)^. 
Furthermore, the canonical metric on dP^ is simply the one induced from the Fubini-Study 
metric on (CP^)^. 

While the action of the symmetries in symplectic coordinates is geometrically clear, we 
want to make explicit their description in complex coordinates. Note that the action of R2 
on the complex coordinates Ui is given by the inverse transpose of R2 since u = dg/dx. Thus 
R2 acts as (^1,^2) I— ^ iu2,U2 — ui). Similarly, Ri which is its own inverse transpose, acts 
by {ui,U2) ^ {u2,ui). More explicitly, the relation between the complex and symplectic 
coordinates for the canonical potentials is given by 



1 + Xi 1 + Xi + X2 
1 — Xl 1 — Xl — X2 



(3.51) 



1 + X2 1 + Xi + X2 
1 - X2 1 - Xi - X2 



To invert these relations involves solving a cubic equation in one of the Xj. 

There are six afhne coordinate patches associated to the cones formed by pairs of neigh- 
boring rays Va- We denote the cone formed by the ray Va and Va+i to be aa and the coordinate 
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system on a a by (^ai Va)- The complex coordinates on the patches a = 1, . . . , 6 are 



(6,»?2) 



(y,xy-i) 
(x-i,y) 

(x,y-i) 



(3.52) 



Note that rja = 1/^a+ii and that when = 7?a+ii then ^a+i = rja- Because of these relations, 
we can consider an atlas on dP^, where for each a a, we restrict ^ 1 and ^ 1- These six 
polydisks Pa tile dP^. 

In the symplectic coordinate system (xi,X2), our atlas divides the hexagon up into six 
pieces. The boundaries are given by the conditions e"^ = 1, e"^ = 1, and e"^~"^ = 1 or in 
symplectic coordinates, the three lines 2xi + X2 = 0, xi + 2x2 = 0, and xi — X2 = 0. 

The action of Dq maps one Pa into another. For example in patch 3 since X and Y are the 
exponentials of ui and U2-, R2 acts on our coordinates by sending {X^"^, Y) (y~^, X~^y), 
mapping patch 3 into patch 4. 



Figure 2: The polytope P for dP^. The dashed lines correspond to the edges of the unit polydisks: 
2.T1 + X2 — xi + 2x2 ~ 0, and xi — X2- The red curved lines correspond to setting and r/i equal 
to two instead of one in the appropriate coordinate systems. 

4. Ricci flow 
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where we have set A = 1. Note that, smce the Ricci form is a closed two-form, a Kahler 
metric remains Kahler along the flow. Furthermore, if [uj] = [TZ] for the initial metric, then 
the flow will stay within that Kahler class. If the flow converges, the limiting metric will 
clearly be Einstein. A recent result of Tian and Zhu [14] implies that, on dP^, starting from 
any initial Kahler metric obeying [uj] = [TZ], the flow will indeed converge to the Kahler- 
Einstein metric. Numerical simulation of the flow can therefore be used as an algorithm for 
finding the Kahler-Einstein metric. 



Thanks to the toric symmetry, (4.1) can be written as a parabolic partial differential 
equation for a single function in 2 space and 1 time dimensions, and therefore represents a 
simple problem in numerical analysis. Specifically, according to (3.22), we can write ( |4.1| ) in 
terms of the Kahler potential f{u): 

Of 

^ = lndetF,j+2/ + c. (4.2) 

Here c, which is independent of u but may depend on t, can be chosen arbitrarily. In particular, 
the zero mode of / (which is pure gauge) is clearly unstable according to ( f4.2D , and c is useful 
for controlling that instability. In principle we could also add a term 7 • n to the right-hand 



side. As explained in Section 3.8, however, due to the symmetry of the polytope for dP^, 



we know that 7 = in the solution to the Monge- Ampere Eq. (3.24). In the symplectic 
coordinates, ( [4. 21) takes the form 



|.lnde.<.,..(-..| 



Indet Gij + 2{-x - ^+g] -c. (4.3) 



Note that, because the Kahler potential is t-dependent, the mapping relating u to x is as well. 
To derive ([4.3|), we used the fact that, since / and g are related by a Legendre transform, 
dg/dt\x = —df/dt\u- Note also that ([4.3|) does not quite describe Ricci flow in the symplectic 
coordinate system; rather it describes Ricci flow supplemented with a t-dependent diffeomor- 
phism (described by the flow equation dg^^/dt = —2R^y + 2g^y + 2V(^^j^-)). Under pure Ricci 
flow, symplectic coordinates do not stay symplectic, since the Ricci tensor ( 3.3C| ) is not the 



Hessian of a function, hence the necessity of supplementing the flow with a diffeomorphism. 

The two flows ( |4.2| ) and ( |4.3| ) were simulated by two independent computer programs, 
which yielded consistent results. In both cases we represented the (Kahler or symplectic) 
potential using standard real space finite differencing, and therefore obtain the resulting 
approximation to the geometry in the form of the potential at an array of points. (In the 
next section we discuss how to present this information more compactly using polynomial 
approximations.) We relegate technical details of these implementations to Appendix ^. 
When using finite difference methods it is important to show that one obtains a suitable 
convergence to a continuum limit upon refinement of the discrete equations, and we present 
tests that confirm this, and estimate errors in the same appendix. 
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4.1 Complex coordinates 

We now discuss in more detail the implementation of Ricci flow using complex coordinates. 
Our canonical Kahler potential /can is defined in P°, the interior of the hexagon. To work 
on our atlas ( 3.52| ), we use Kahler transformations to modify the Kahler potential so that / 
is smooth on each coordinate patch, including the two relevant edges of the hexagon. This 
Kahler transformation must also respect the ?7(1) isometries so must only depend on the 
Ui. The only such Kahler transformations are affine linear combinations of the Uj. (For a 
function h{ui+i9\,U2+i02) depending only on the Ui, the condition on Kahler transformations 
ddh = reduces to g^^g^ =0.) In the interior of the hexagon, we can thus modify / without 
modifying the metric by adding a linear combination of ui and U2- In symplectic coordinates 



1 ^ 

■ K = - ^^Va - K \n{Va ■ X + I) 



a=l 



where k G M^. 

We find that for any Kahler potential /, regardless of coordinate patch. 



detK 



(-2 2 



We can rewrite ( [4.2[) , yielding 



^ J \r] dr] drf 



idf , d^f 



(4.4) 



dt 



^ 9^ J \r] dr] drp' 



\d^drj J 



+ 2(/ + lnC + lnj?) +c . (4.5) 



This shift of / by logarithms is a Kahler transformation. In each coordinate patch, we can 
write In^a + ^ri-Va = u ■ Ka where Ki = (1,0), K2 = (0,1), K3 = (-1,1), K4 = (-1,0), 
K5 = (0,-1), and kq = (1,-1). Moreover, for these choices of Ka, fa = f + u-Kais well 
behaved everywhere inside Pa including the edges. 

We simulate Ricci flow using a Kahler-Einstein potential fa defined on a neighborhood 
Ua of Pa, Ua = {{^a,Va) '■ ia,Va < L, L > 1} which Satisfies (^). As an initial condition, 
we take fa{t = 0) = /can + In^a + Inr/a- Given two patches, a and 5, then for a point in the 
overlap, p G UaCiUb, 

fb{p) = fa{p) + U- {Kb- Ka) . (4.6) 

These Kahler transformations are consistent with the canonical potential and thus fix the 
same Kahler class. There is then an element R £ Dq such that R{Pb) = Pa which relates fa 
and fb- In particular, if p in is in Ub and Rp is in Ua- 



fbip) = faiRp) . (4.7) 

Thus, with these quasiperiodic boundary conditions along the interior edges ^a = L and 
Tja = L, we can determine / by working solely on the patch Ua- 
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In addition to quasiperiodic boundary conditions along the interior edges of Ua, close to 
the exterior boundary of [/q, we use the condition that the normal derivative to the boundary 
must vanish. These Neumann boundary conditions arise because in complex coordinates, 
approaching the exterior boundary should be like approaching the center of the complex 
plane. 

The Ricci flow in this domain was represented using second order accurate finite differ- 
encing with various resolutions up to a grid J-jj of 250 x 250 points. This discretization, 
its convergence to the continuum and error estimates at this resolution are discussed in the 
Appendix]^, and data presented in this section is given either from extrapolating to the con- 
tinuum or using this highest resolution. In particular we estimate that the Kahler potential 
computed at this resolution is accurate pointwise to one part in 10^. 

We noted earlier that while Ricci flow may converge, the zero mode of the potential / is 
not guaranteed to converge since it does not appear in the metric. For convenience we wish 
to obtain a flow of / that does converge, so we promote the constant c to be time dependent 
along the flow, but still a constant on the geometry (i.e. we introduce a flow dependent 
constant Kahler transformation), and choose c such that the value of J^n/2,n/2i where N is 
the grid size, does not change. At large times, c tends to a constant. At the end of the flow, 
we found it convenient to set c = by adding an appropriate constant value to the grid J-. 

4.2 Symplectic coordinates 

We move on to discuss the implementation of Ricci flow based on symplectic coordinates. As 
discussed at the beginning of the section, although we work in symplectic coordinates, we use 
Ricci flow defined by the complex coordinates; since the relation between the two coordinate 
systems moves around as the metric changes, in the symplectic coordinates this is Ricci flow 
plus diffeomorphism. 

In order to deal with the boundary conditions on the edge of the polytope, it is useful 
to work with h rather than g (recall that g = gc^n + h and Qcan is the canonical symplectic 
potential). In terms of h, ( [4. 3D becomes 

f)h / rP'h \ f r)h \ 

m = (!" + + ft? + ^) + - l^'^' 

where 



In det + 2 ( -X • %^ + 5can ) = ln(Li + 1^ + 1^), (4.9) 



a 



dx 

2 



1 + zt 



-1 1 + zf 



(4.10) 



Li = hh = l-xi, L2 = hh = l-xi, Li = l2h = l-{xi+X2Y. (4.11) 

Due to the hexagon's Dq symmetry, it is sufficient to simulate the flow within the square 
domain < xi < 1, — 1 < X2 < 0. Ricci flow in this domain was represented using second 
order accurate finite differencing with various resolutions up to a grid of 400 x 400 points. 
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The differencing, continuum convergence and errors at this resolution are discussed in the 
Appendix ^ with data presented in this section given either from extrapolating to the con- 
tinuum or using this highest resolution. In that appendix we estimate that the symplectic 
potential computed at this resolution is accurate to one part in 10^ at a given point. 

As for the complex coordinates, the zero mode of the symplectic potential is pure gauge 
and does not converge. The constant c was promoted to depend on flow time and chosen to 
keep h{0, 0, t) = 0. 

4.3 Results 

As predicted by the Tian-Zhu theorem, the metric converges smoothly and uneventfully to 
the Kahler-Einstein one. In the symplectic implementation, the Ricci flow was simulated 
starting with a variety of initial functions ho{x) = h{x,t = 0) (always corresponding, of 
course, to positive-definite initial metrics). Three examples are shown in Fig. ^. For every 
initial function investigated, the flow converged to the same fixed point h^ix) = h{x,t = oo), 
which necessarily represents the Kahler-Einstein metric. The exponential approach to the 
fixed point is controlled by the scalar Laplacian; this will be discussed in detail in the next 
subsection. 

The final complex and symplectic potentials found by the independent implementations 
are plotted in Figs. ^ and ^, along with the respective canonical potentials. The results are 
plotted in the fundamental domain actually simulated, and one should use the Dq symmetry 
to picture the potentials extended over the whole domain. To compare the Kahler potential 
fE{u) computed in complex coordinates with the symplectic potential gE^x) computed in 
symplectic coordinates, we numerically performed a Legendre transform to obtain a Kahler 
potential f^{u) from gE{x). Plotting — /e in Fig. ^, it can be seen that they differ by less 
than 5 x 10"^. (In fact the agreement may be slightly better as there is likely some error 
introduced in doing the numerical Legendre transform.) Thus our results appear to agree to 
around the same order that we believe they are accurate. 

In order to get some feeling for the form of the Kahler-Einstein metric, and how it 
compares to the canonical one, it is helpful to plot some curvature invariants. Of course, any 
invariant depending solely on the Ricci tensor, such as the Ricci scalar, will be trivial, so we 
need to go to invariants constructed from the Riemann tensor. (Expressions for the Riemann 



tensor are given in Section |3.4| above.) For example, the sectional curvature of the xi-X2 plane 
(at fixed 9i), which is -Ra;iX2a;iX2/ det(Gjj), is plotted for the canonical and Einstein metrics 
in Fig. 0. Also of interest is the Euler density 

e = ^ (i?' - ^R^.uR^" + R^.upxR'"''^) , (4.12) 

which integrates to the Euler character of the manifold, which is 6 for dP^. The first two 
terms inside the parentheses cancel in the case of an Einstein metric. Fig. ^ shows 47r^e for 
the canonical and Einstein metrics. The factor of 47r^ takes account of the coordinate volume 
of the fiber. Recalling that = 1 in symplectic coordinates, the plotted quantity should 
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Figure 3: The Ricci flow using symplectic coordinates starting with (top) the canonical metric 
{ho — 0), (middle) ho = —^{^i + + X1X2), (bottom) ho = (x^ + X2 + xiX2)^. Left side: g ^ ^can + h 
versus xi along the line X2 — 0, plotted at intervals of O.f units of Ricci flow time (t = is the top 
curve). Right side: g versus t at the (arbitarily chosen) point {xi,X2) — (0.3,-0.2). We see that all 
three initial conditions converge to the same fixed point. 

integrate to 2 over the plotted region, which covers one-third of the polytope. This can easily 
be checked in both cases by numerical integration. 

4.4 Laplacian eigenvalues 

An important geometric quantity is the spectrum of the scalar Laplacian. Here we illustrate 
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Figure 4: a) The Einstein Kahler potential /e as a function of ^ and r/; b) the canonical potential 
/can($,'7); and c) /E(^,r?) - /can(^,r?)- 

a simple method to compute low-lying eigenfunctions. The natural flow associated with the 
scalar Laplacian is diffusion, and the late time asymptotic behavior of the diffusion flow 
is dominated by the eigenfunction with lowest eigenvalue. Hence simulating diffusive flow 
on the dP'i geometry found, and extracting the asymptotics of this flow allows the lowest 
eigenfunction to be studied. We may classify the eigenfunctions under action of the Dq and 
isometries. Since the flow equation is invariant under these symmetries, if we start 
with initial data that transforms in a particular representation, the function at any later time 
in the flow will remain in this representation. For simplicity we will focus on eigenfunctions 
which transform trivially, but obviously the method straightforwardly generalizes to compute 
the low-lying eigenfunctions in other sectors. As for the Ricci flow, the flow does not depend 
on second normal derivatives of ip at the boundaries of the hexagon domain, and hence we 
do not require boundary conditions for if) there, except to require it to be smooth. 

The lowest eigenfunction of — Ae in the symmetry sector we study is -0 = constant 
which has zero eigenvalue. We are interested in the next lowest mode which has positive 
eigenvalue and non-trivial eigenfunction, denoted ipi{x) with eigenvalue Ai. Then we consider 
the diffusion flow on our del Pezzo solution, 

= AEV(t,^) (4.13) 

and start with initial data for tjj that is symmetric and will hence remain symmetric. At late 
times, the flow will generically behave as, 

0(t, x)=i)Q + Vi(^)e"^i* + 0(e^^2*) (4.14) 

where V'o is a constant, corresponding to the trivial zero eigenmode, and A2 is the next lowest 
eigenvalue A2 > Ai. Waiting long enough and subtracting out the trivial constant, the late 
flow is given by ipi, the eigenfunction we wish to compute. 

In Fig. P we plot the log of 0, 0) — ipQ as function of the flow time t for 3 different 
initial data. Once the higher eigenmodes have decayed away, we clearly see the flows tend to 
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Figure 5: Top left: The canonical symplectic potential gcan- Top right: The Einstein symplectic 
potential = ffcan + ^e- Bottom: h^- These are plotted in the range 0<a;i<l,— l<a;2<0, 
which is one-third of the hexagon; the values on the rest of the hexagon are determined from these by 
its Dq symmetry. 



the same exponential behaviour. We estimate this eigenvalue by fitting the exponential decay 
as Ai = 6.32. In Fig. 10 we plot the eigenfunction 'ipi{x), normalized so that 1^1(0, 0) = 1. 
Note that, as expected, for different initial data we consistently obtain the same function. 

This lowest eigenvalue and eigenfunction can also be obtained from the approach to 
the fixed point of the Ricci flow. For concreteness let us work in symplectic coordinates; 
corresponding expressions will hold in complex coordinates. Expanding h about its fixed- 
point value, 

h = hE + Sh, (4.15) 
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1 



Figure 6: The difference between the final Kahler potential computed in complex coordinates and 
the Legendre transform of the final symplectic potential computed in symplectic coordinates. 




Figure 7: The sectional curvature of the X1-X2 plane at fixed angle 9i, for the canonical (left) and 
Einstein (right) metric. The two are quantitatively different but qualitatively similar. 



the flow equation ( [4. 8]) becomes, to first order in 6h 

d6h 



(Ae + 2)5h - 5c, (4.16) 

where 6c depends on how c is chosen. Since we used initial conditions for the Ricci flow that 
respected the Dq symmetry, we should find that the potentials approach their fixed point 
values the same way as ip above, with a shift of 2 in the exponent: 

/i(t,x) =/iE(x) + (Vi(x)-^i(0))e-(^i-2)* + ... ; (4.17) 
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Figure 8: 47r^e, where e is the Euler density, for the canonical (left) and Einstein (right) metrics. 
The Euler density is distributed differently but integrates to the same value — 2 over this region, and 
6 over the whole polytope — for any two metrics. It is more evenly distributed for the Einstein metric 
than for the canonical metric, but still peaked at the vertices of the polytope. 

Log[(A(t,0,0)-iAo] 




Figure 9: Decay of tp under diffusion towards a constant as a function of diffusion time t. The log 
of ^l;{t, 0,0) — ^po is plotted, the slope giving the eigenvalue for the lowest (non-constant) symmetric 
eigenf unction. The 3 curves correspond to 3 different initial profiles for tp although we see that the 
decay quickly becomes dominated by the lowest eigenmode. 

the eigenfunction is shifted by a constant because c was chosen to keep h{t, 0) = along the 
flow. Our numerical flows confirm this expectation. The corrections in ( [4.17| ) involve both 
higher eigenvalues of the Laplacian and higher-order effects in 6h. 
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1-1 



Figure 10: The lowest eigenfunction of the scalar Laplacian which transforms trivially under the Dq 
and J7(l)^ symmetries. 



5. Symplectic polynomials 



In the Ricci flow simulation in symplectic coordinates discussed in Section 4.2, the function 
h{x) — which encodes the symplectic potential and therefore the metric — was represented 
by its values on a lattice of points in xi, 0:2. In this section we will discuss a different way to 
represent the same function, namely as a polynomial in xi,X2- Since the solution to the 
Monge- Ampere equation is a smooth function, it can be represented to good accuracy with 
a vastly smaller amount of data in this way: a few polynomial coefficients as compared to 
values on thousands of lattice points. Furthermore, quite independent of the solutions found 
in the previous section, the problem of finding an approximate solution to the Monge- Ampere 
equation can be expressed as an optimization problem for the polynomial coefficients; we will 
use this fact to develop a third algorithm in the following section that is quite different in 
character from Ricci flow. 

It is interesting to note that the metrics obtained from polynomial expressions for h{x) 
are the symplectic analogues of the so-called "algebraic" metrics on Calabi-Yau manifolds 
that have been used for numerical work by Donaldson [3] and Douglas et al. [4,5]. The 
algebraic metrics, which are defined for a Calabi-Yau embedded in a projective space, have a 
Kahler potential that differs from the induced Fubini-Study one by (the logarithm of) a finite 
linear combination of a certain basis of functions, namely the pull-backs of the Laplacian 
eigenfunctions on the embedding projective space. This is a generalization of the usual 
strategy of representing a function by expanding it in a basis of Laplacian eigenfunctions 
(such as Fourier modes); since the eigenfunctions on the Calabi-Yau depend on the metric 
that one is trying to find, one instead uses the eigenfunctions on the embedding space, which 
are known in closed form (and are indeed very simple). In our case, we consider the embedding 



of dP^ in (CP^)3, which, as discussed in Section |3^ , is described in symplectic coordinates by 
the equation xi + X2 + = (where the are the symplectic coordinates on the respective 
CP^ factors). The first n eigenspaces of the Laplacian on (CP"'^)^ (with respect to the Fubini- 
Study metric), restricted functions that are invariant under the U{1)^ isometry group, are 
spanned by the monomials in xi,X2,X3 up to order n. This is precisely the basis of functions 
we use to expand the difference h between the symplectic potential g and the induced one 

5can- 

We now describe some fits to the numerical solutions of the last section in terms of 
polynomials up to sixth order in xi and X2, and quantify how well those polynomials do in 
solving the Einstein equation. These small polynomials likely provide sufficiently accurate 
approximations to the Einstein symplectic potential for most purposes, while at the same 
time being more tractable than the full numerical data for analytic calculations. 

We begin by noting that, since He is invariant under the hexagon's Dq symmetry group, 
it is sufficient to consider invariant polynomials. As shown in Appendix every invariant 
polynomial can be expressed in terms of the two basic invariant polynomials. 



U = xl + X1X2 + xl , 



V = XiX2{xi + X2) 



(5.1) 



We simply do a least-squares fit of the polynomial coefficients to the lattice values of He 
obtained in the Ricci flow in symplectic coordinates, that is, we minimize 



1 



a 



VdPs JdP'i 



\^{Kt - hEY 



1 



VdPa JdPs 



Vg {hfit - /ie) 



(5.2) 



(Any constant difference between /igt and He is irrelevant). At successive orders in x we find 
the following fits: 



^fit 


a 


P 









0.06 


0.5 


-0.24C/ 






IQ-^ 


0.1 


-0.2214f7 - 


0.0215f72 




10-4 


0.03 


-0.22412^7 


- 0.01450C/2 


- 0.00521C/3 + 0.00734y 


10-5 


0.007 



(5.3) 



In each case we have written only the significant digits of the coefficients. ^ Independently of 
our numerical result He, it is useful to know how far the metric corresponding to hfn deviates 



from being Einstein. In Fig. 11, the pointwise rms deviation of the eigenvalues of the Ricci 
tensor from 1, 



(5.4) 



These digits do not change between the run with 200 lattice points and the run with 400 lattice points 
(except the last digit of each coefficient in the sixth-order approximation), and are therefore presumably equal 
to their continuum values. 



- 26 - 



is plotted for these four functions. The global rms deviation from being Einstein, /3, where 



1 



(5.5) 



KiPs JdP'i 

is also shown in the table above. As expected, each successive order gives a substantially 
better approximation, and a metric that is substantially closer to being Einstein. 




Figure 11: D = \J \{Rfiu — .9p;y)^ versus Xi,X2 for the four polynomial functions listed in ( |5.3| ). The 
deviation from being Einstein is most significant along the edges of the hexagon in the first two cases, 
and at the corners in the second two. 

For eigenfunctions of the Laplacian, we may perform the same invariant polynomial fits 
as we did for h; at quadratic order we find ^/^i ^ 0.985 — 2.37U. 

We regard table ^.3| as a key result of this paper. The last line of the table provides in an 
extremely compact, analytic form, an approximation to the true symplectic potential on dP^. 
It deviates from the true potential pointwise by at most ~ 0.1%, and satisfies the Einstein 
condition well within the hexagon, giving at most a 10% error near the hexagon corners as 
measured by the pointwise rms deviation of the Ricci tensor eigenvalues, defined above. 
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6. Constrained optimization 



In the previous section we used the results of Ricci flow to find a polynomial approximation 
to h(x), the smooth part of the Kahler-Einstein symplectic potential (recall that h = g — 
Qcan)- Here we instead search for a polynomial approximation to h using the Monge- Ampere 
equation directly. A simple approach would be to consider the space of polynomials of some 
given order, and minimize an error function built from the Monge-Ampere equation on that 
space. However, if one wishes to obtain high accuracies, one needs to go to high orders, 
and then this brute force approach rapidly becomes intractable due to the large number 
of polynomial coefficients and the difficulty of searching in a high dimensional space. We 
therefore take advantage of the analytic properties of the Monge-Ampere equation to constrain 
the polynomial coefficients, by requiring the polynomial to solve it order by order in Xj. As we 
will see, this leaves only a small number of undetermined parameters, dramatically simplifying 
the error function minimization. We now explain the details of the method. 

We begin by noting that the exact solution h^ix) to the Monge-Ampere equation is an 
analytic function of the Xi, which can be seen in a couple of ways. As noted, the symplectic 
coordinates are eigenfunctions of the Laplacian. Hence, in these coordinates the Ricci curva- 
ture operator has the same character as it does in harmonic cooordinates: it is actually an 
elliptic operator. Because the Einstein equations are analytic in the metric, they will have 
analytic solutions. Somewhat more directly, the Monge-Ampere equation we are solving is 
elliptic at the Kahler-Einstein potential and analytic, and so the solution will be analytic. 

Before we constrain the polynomial approximation to h{x) using the Monge-Ampere equa- 
tion, we constrain it by imposing the hexagon's Dq discrete symmetry group. As discussed 
in Section ||, any invariant polynomial in Xi can be written in the form 

where U and V are given in Eq. ( |5.1| ). To eighteenth order in xi and X2 we write the series 
as follows 

h = A0 + A1U + + + ^4C/4 + A^U^ + AqU^ + ArU^ + AsU^ + AgU^ + ... 
+ ViBo + BiU + B2U^ + B^U^ + BiU^ + B^U^ + BqU^ + . . .) 
+ (Co + CiU + C2U'^ + CsU^ + ...) + V^ (Do + ...) + ... . (6.2) 



Plugging this series into the Monge-Ampere equation ( 3.27 ) (with 7 = and A = 1), 
yields constraints on the Cij that relate the Cij with j > to the q^q- To make the expressions 
a little simpler, we introduce a new constant a: 

= -lln3-lna . (6.3) 

We worked out the relations up to order 18 in xi and X2- The first relation is that Ai = —Izba. 
The numerical Ricci flow results are consistent only with the plus sign. The next few relations 
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The power of this method is that since the Monge- Ampere equation determines many of 
the coefficients, when we determine the remaining ones by minimizing an error function, then 
at a given order there are far fewer parameters to solve for. Truncating at eighteenth order 
in xi and X2, we need fit only 4 parameters Aq, Bq, Cq, and Dq to find an approximation to 
h. In contrast, minimizing an error function using the most general unconstrained eighteenth 
order polynomial given in equation |6.2| involves searching a 22 dimensional space. 

In principle if we took an arbitrarily high order expansion, then about the origin of the 
hexagon, in the region where the series for h converges, we would solve the Monge-Ampere 
equation precisely. However we see that we still have undetermined constants in the series, 
and these correspond to the fact that we must provide boundary conditions to determine a 
solution fully. 

A posteriori, our numerics strongly suggest that h converges everywhere in the interior 
of the Delzant polytope. Constraining h to have the correct behavior on the boundary of the 
polytope should completely determine the remaining constants in the power series expansion. 
Hence we fit the remaining parameters by requiring Eq. ( |3.27D be satisfied at the boundary of 
the hexagon. As emphasized in the discussion around Eq. ( 3.29| ), the boundary conditions do 
not need to be supplied separately — they are enforced by Eq. ( |3.25D itself. For dP^, along 
xi = 1, Eq. ( |3.25 ) reduces to the boundary condition 



= 1 - 2x2 - 2x1 + (1 - 4)i'^ - (1 + X2f))hx2x, - exp [2(/i^, + X2h^, - h)] = 



(6.4) 

To find this expression, we have assumed that h is smooth at the boundary xi = 1. At 
the corners X2 = and X2 = —1, Eq. ( |6.4D reduces further to h = hxi and h = —hx2 + 
respectively. Explicitly, we determine the remaining parameters (^o through Do in the 
eighteenth order truncation) by minimizing 



M 



p 



\BC{p)\' 



summed over twenty equally spaced points along the boundary xi 
more terms in the series, 



(6.5) 

1. We find, as we include 



order 


M 


a 


Bo 


Co 


Do 


2 


0.02 


0.757 








6 


9 X 10-^ 


0.7753 


0.011 






10 


3 X 10-^ 


0.77616 


0.00508 






14 


5 X 10"*^ 


0.776226 


0.00480 


-0.00055 




18 


2 X 10"^ 


0.776235 


0.004781 


-0.00015 


0.004 



(6.6) 
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The full expression for the 18th-order polynomial is given in a Mathematica notebook available 
for download at the websites [21]. In Appendix]^ we show the dependence of h at various 
locations in the polytope as a function of the number of terms taken in the expansion. We 
see that convergence for h is fast — apparently faster than polynomial — in the number of 
terms, and in particular for everywhere tested within the hexagon, and also on its boundary, 
we see convergence. In particular we see no sign of poor behaviour near or on the boundaries 
of the hexagon. We also see that the values the series converges to are in excellent agreement 
with the continuum extrapolated values of h found from the Ricci flow method and detailed 
in Appendix From these data we estimate that the potential given by the eighteenth 
order expansion differs from the true solution by approximately one part in 10^, and hence 
is comparable in this respect to the 400 x 400 Ricci flow result. 

The figure of merit M does not give a very good indication of the degree of accuracy of 
our fit globally. To understand how well we are doing globally, we use the same local estimate 
of error as in the previous section, D = ^J \{R^i^ — g^v)"^- We find that the maximum value D 
attains in the domain decreases with each increase in order of the expansion. The maximum 
is found on the lines connecting the origin to the hexagon vertices, and hence in Fig. |l^, we 
plot this error estimate along one of these lines, L'(xi,0), for the 6th, 12th and 18th order 
polynomial expansions. As expected the error is smallest at the origin, and most error is 
localized near the boundaries. Since the error is rather localized near the hexagon vertex, we 
have plotted this error against In (1 — xi) to demonstrate that it is indeed finite at the vertex. 
We see that while the symplectic potential taken pointwise may be accurate to 1 part in 10^ 
as stated above, since the error is localized in the hexagon corners, the quality of the solution 
is a little worse in these regions. We see for the 18th order approximation that the error in 
the Einstein condition, estimated by D, is about one part in 10'^ at the vertex. 
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Figure 12: D, the error in the Einstein condition, along the line X2 = from the origin xi = to the 
hexagon corner xi = I where the maximum of D in the hexagon domain occurs. The left plot is for 
a 6th order expansion, the middle is 12th order, and the right is 18th order. We see the error in the 
Einstein condition is quite localized near the hexagon corner (hence we plot against In (1 — xi)), but 
remains finite there, and decreases everywhere with increasing numbers of terms in the expansion. 
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6.1 Laplacian eigenvalues 



We can use this 18tli order fit to extract the lowest eigenvalues (of eigenfunctions invariant 
under the Dq action) on our manifold. Note that in coordinates where the metric is analytic, 
eigenfunctions of the Laplacian are analytic: elliptic equations with analytic coefficients have 
analytic solutions. Thus we can play a very similar game, expressing the eigenvectors as a 
local power series near the origin of the hexagon in the U and V variables: 



i^{U, V) = ^ + XiU + + XsU^ + YiV + ... 



(6.7) 



We have chosen to normalize ^(0,0) = 0.1. We could in principle have used the differential 
equation to constrain some of the Xi and Yi, but we did not, aiming for a fit whose errors 
are more evenly distributed over the hexagon. We fit a hundred equally spaced points in a 
square domain < xi < 0.9 and —0.9 < X2 < and minimize 



(6.8) 



as a function of A and the X^ and Yi. By searching for successive local minima of M^, we 
can extract successively higher eigenvalues. We find 



order 




Xi 


X2 


^3 




Ai 


2 


0.006 


-0.239 








6.27 


4 


0.002 


-0.246 


0.011 






6.325 


6 


3 X 10"^ 


-0.245 


0.006 


0.006 


-0.024 


6.322 



(6.9) 



order 




Xi 


X2 


^3 


Yi 


A2 


4 


0.5 


-0.69 


0.79 






17.4 


6 


0.008 


-0.67 


0.70 


0.16 


-1.28 


17.2 



(6.10) 



For the lower eigenvalue 6.32, note that the ratio of the first two coefficients —2.39 is in 
reasonably good agreement with the corresponding ratio —2.37/0.985 = —2.41 determined in 
Section ^. 

The choice to minimize in the domain < \xi\ < 0.9 was a compromise that requires 
some justification. First, since ( ^.25 ) enforces its own boundary conditions, minimizing 
close to the boundary xi = 1 is enforcing the boundary condition to first order in xi — 1. 
Second, the power series approximation for tp experimentally does not appear to have good 
convergence properties near the boundary. Experimentally, the best values for the coefficients 
of the truncated power series (in the sense of agreeing with the coefficients of the power series 
itself) are obtained by making a compromise between minimizing over a set of points that 
extends to the boundary and minimizing over a set of points for which the power series has 
good convergence properties. 
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6.2 Harmonic (l,l)-forms 

As noted in Section harmonic (1, l)-forms and eigenfunctions of the Laplacian must satisfy 
a very similar equation. We end this section with a computation of the harmonic (1, l)-form 
9 a- Unhke the eigenfunctions computed above, 6a does not transform trivially under Dq. 
Thus, we assume that has a more general expansion of the form 



IXa = ln(l + • a;) + ^ CnmXiX^ 



(6.11) 



We use the same least squares approach as above, minimizing 

Me = ^\/\^la 



const 



Because of the explicit xi <-> X2 symmetry, we start with a 
to sixth order in xi and X2, we find 



2 and set d 



(6.12) 
Cmn- Fitting 



ClO = 


-0.2250 










C20 = 


0.0638 


Cll = 


0.0311 






C30 = 


-0.0301 


C21 = 


0.0059 






C40 = 


0.0126 


C31 = 


0.0005 


C22 = 


-0.0073 


C50 = 


-0.0150 


C41 = 


-0.0245 


C32 = 


-0.0240 


C60 = 


0.0088 


C51 = 


0.0223 


C42 = 


0.0281 



(6.13) 



C33 = 0.0196 



The value of Mq ~ 10~^ at the minimum implies an average error of 10^'^ at each of the 100 
points. Note the error gets much worse outside the fitting domain Ixjl > 0.9. The philosophy 
in this section is similar to that in the discussion of eigenfunctions: we are attempting to find 
more accurate values of the Cij rather than attempting to minimize the global error. 
The fit also yields F^Wij = 0.6672, consistent with our expectations. We know that 



1 



UJ = - 



2 ^ 

a=l 



(6.14) 



Clearly F^^Fij = 2. By the dihedral symmetry group, the value of F^'Wij should be independent 
of a. We conclude that 

{eahr^ = -. (6.15) 



From 02, we can reconstruct the other 9a by applying the Dq group action. To test how 
good our approximation to 02 was, we computed = Q^^Qij where Qij = J2a(^a)ij using 
our best fit for 02- Since Qij should be 2Fij, should be approximately eight. A plot of @^ 



is shown in Fig. 13 



For the purposes of the KT solution described in the Introduction, we need a such that 
0^^Fij = 0. From the preceding discussion, any linear combination of the form '^a^a^a such 
that Ca = will have this property. We also require that ^6* = —0. In fact, the condition 
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Figure 13: The value of 9^ for our 6th order fit. 



Xla Ca = enforces anti-self-duality. The reason is that the Hodge star treated as a linear 

operator acting on the space of harmonic (l,l)-forms has signature (H ). We know that 

-koj = LO] thus any (l,l)-form orthogonal to a; must be anti-self-dual. In general, the numerics 
suggest that for such a 0, O^Wij will be a nontrivial function of both Xi and thus that solving 
for h{p) requires solving a PDE in three real variables. 



7. Discussion 



In this paper we have described three different methods to find the Kahler-Einstein metric 
on dP^. All three methods exploit the Kahler and toric structures of the manifold, allowing 
us, using modest computing resources, to compute the metric in both Kahler and symplectic 
coordinates to an accuracy of one part in 10^. The results of the different methods are 
consistent to within that error. We expect that this accuracy is sufficient if one wishes to 
compute geometric quantities for either physical or mathematical applications, and we have 
made available the data along with Mathematica notebooks to allow manipulation of these 
results [21]. 

We noted that, for a lesser accuracy of one part in 10^, a simple expression for the smooth 
part h of the symplectic potential, g = gcan + h, already provides such an approximation, and 
we repeat it here: 

h{xi,X2) = -0.22412C/ - 0.01450[/2 - 0.00521C/3 + 0.00734^, 

U = xf + XiX2 + X^^ V = x\x\{xi + X2f' ■ 



where ^can is given in equation 3.15. The resulting metric satisfies the Einstein condition 



everywhere to better than 10% as discussed in Section ||. 
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Simulation of Ricci flow has proven to be an effective way to solve the Einstein equation. 
We have found that implementing the flow is a little simpler in symplectic coordinates: the 
domain is naturally compact; the symmetries are more manifest; and the boundary conditions 
are simpler. Our codes (which were not optimized for speed) converged in a few hours for 
the highest resolutions. For higher accuracy than attained here, one could optimize the flow 
simulation, for example by taking more advantage of the discrete symmetries than we have 
done. More generally, Ricci flow simulation, using an explicit finite differencing method as we 
have done, can be thought of as a particular iterative scheme for solving the Monge- Ampere 
equation. If one is interested only in solving that equation, and not in accurately simulating 
Ricci flow, then this scheme could be modified to improve speed. For example, by replacing 
the Jacobi-type updating method by a Gauss-Seidel method, one obtains a faster algorithm 
(experimentally, 50% faster in complex coordinates). To obtain a parametric improvement in 
speed would likely require a non-local modification such as multi-grid. 

The constrained optimization approach we have demonstrated uses the symplectic poly- 
nomials, reducing the size of the search space by solving the Monge- Ampere equation order 
by order in Xi, and has proven very powerful. It is as accurate as the Ricci fiow results, 
but is quicker. One drawback is that the hexagon origin is singled out as the point where 
the solution is best, and the error in the solution becomes tightly localized at the corners 
of the hexagon. Most computational time is invested in determining the constraints in the 
series expansion, and this algebraic problem gets worse the more terms that are included in 
the expansion. However, once one has this solution, the numerical minimization of the error 
function is simple. 

The two methods are complementary in the sense that for Ricci flow the time is spent in 
numerically computing the flow, whereas for the optimization the time is spent algebraically 
computing the expansion of the potential. The principle advantage of Ricci flow is that the 
method is very general, and while it benefits from the Kahler and toric structures it certainly 
applies to more general problems which do not possess them. It is not clear how widely 
applicable the constrained optimization approach is, as it likely works due to the special 
properties resulting from those mathematical structures. However, it would be interesting to 
investigate its application to other situations. It would also be interesting to compare these 
approaches, particularly the constrained optimization, with Donaldson's method [3]. 
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A. Canonical metric on CP^ 

Consider CP^ which has the fan, unique up to S'L(2,Z) transformations, vi = (1,0), V2 = 
(0, 1), and vs = (—1,-1). We choose all the = 1 in gcan to assure that the class of the 
resulting Kahler form is proportional to the first Chern class, 

ffcan = ^ (1 + ln(l + xi) + (1 + X2) ln(l + X2) + (1 - Xi - X2) ln(l - xi - X2) . (A.l) 
From gcan, we can reconstruct the complex coordinates 

i- + xi = 7: 7: — ; 1 + X2 = n — ; (A. 2) 

1 _|_ g2«i _|_ g2tJ2 ' ^ ]_ _|_ g2«i _j_ g2«2 ' ^ ' 

. 3 

1 - Xi - X2 



1 + e2"i + e2"2 
The canonical Kahler potential is then 

/can = I ln(l + e2«i + e^"^) - ^ ln3 - m - tX2 . (A.3) 

The expression gcan satisfies ( 3.25| ) provided c = ln(4/3) and 7 = 0. For CP^, the canoni- 
cal metric is the Fubini-Study metric. In terms of the traditional homogenous coordinates 
{Xi, X2, X-^) on CP^, the Kahler potential is traditionally written, with a different choice of 
normalization of the volume, in the patch X3 7^ 

ln(l + |Xi/X3|2 + \X2/X3\^) . (A.4) 

Thus, we identify IX1/X3I = exp(ni) and IX2/X3I = exp(u2). The canonical metric is always 
Kahler-Einstein for Cartesian products of projective spaces [37]. 



B. Numerical Ricci flow implementation and error estimates 

In this appendix we give technical details of our two finite difference implementations of Ricci 
flow, and also discuss the continuum convergence and estimate errors for the resolutions used. 

B.l Implementation in complex coordinates 

We work on a square domain slightly bigger than a unit polydisk, < < L and < tj < L 
(L > 1) where / is assumed to have Neumann boundary conditions along ^ = and r] = 0. 
Along the internal boundaries, ^ = L and r] = L, we employ a kind of periodic boundary 
condition enforced by the Zg symmetry of the hexagon. We take 



/(e,L) = /(l/L,eL)+lnL, 



(B.l) 



to map the boundary 7] = L for < < 1/L back into our square domain. For 1/L < ^ < 1, 
we take a composition of the above map: 

/(e,L) = /(l/(La,0+ln^'C. (B.2) 

For the £, = L boundary, we take the preceding rules with ry and ^ switched. For the points 
along the boundary with ^ > 1 or r/ > 1, we use the rule 

/(e,r/) = /(l/C,l/??)+21nC + 21nr? . (B.3) 



Note that as indicated in Section [4.1| , due to the Zg symmetry, for any point p = (^,r/), with 
^ > 1 or > 1, the value of f{p) should be related by a Kahler transformation to the value 
of / at a point inside the unit polydisk. 

Although we are only enforcing the Zg symmetry along the boundary of our coordinate 
patch, the symmetry will hold globally. First, our initial potential /can respects the symmetry. 
Second, Ricci flow preserves the symmetry. 

To discretize ([4. 2D, we approximated f{C,v) by its values J^ij on an NxN grid with lattice 
spacing A = L/{N — 1) and used standard second order finite differencing for derivatives. In 
evaluating the first order derivatives at the boundary, we took 



Idl 



^=0 de 



and -— — 
rj orj 



(B.4) 

r)=0 



To impose a discrete version of the boundary conditions, we added extra rows and columns 
along the grid. To impose Neumann boundary conditions, we imposed that J^qj = J^2,j 
and J^jfi = J^i^2- To impose the periodic boundary conditions along the {I,N) and {N,J) 
boundaries, we mapped the the point (I, A'' + 1) or (iV + 1, J) back inside the grid using the 
symmetries and used a bicubic interpolation to compute a best value for /. 

B.2 Implementation in symplectic coordinates 

We represented h on a uniform square lattice in the variables xi, X2- Such a lattice has 
several advantages, aside from simplicity. First, according to ( 3.10| ), these points are also 



spread uniformly according to the measure of any symplectic metric. Second, the lattice is 
itself invariant under the symmetries of the hexagon (since these are elements of GL(2,Z)). 
In fact, in view of this symmetry group, the lattice should in a sense be considered triangular, 
with "edges" running not just horizontally and vertically, but also along the diagonals with 
slope —1. In other words, a given lattice point (xi, X2) has six nearest neighbors: (xi ± A, X2), 
(xi, X2±A), and (xi±A, X2=F'^) (where A is the lattice spacing). First and second derivatives 
were calculated using these nearest neighbor points in a way that was accurate to second 
order in the lattice spacing and respected the hexagon symmetries. Generally speaking, the 
lattice spacing was chosen to be one over an integer, so that the polytope boundaries passed 
through lattice points. In view of the free boundary conditions for h, on these boundaries the 
necessary derivatives were computed by extrapolation, using next-to-nearest neighbor points 
to give third order accuracy. 
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B.3 Simulation of Ricci flow 



The Ricci flow was simulated by an explicit method, with first-order accurate time derivatives 
and a time step of oc A^. The constant of proportionality is of order one, but depends on the 
initial metric since the equation is non-linear. For the canonical choices of initial potentials 
in the complex coordinate case we required a time step ^A^, while in the symplectic case we 
required ^A^. However modifying the initial potential may require a smaller initial timestep. 

The diffusive nature of the flow requires the time step to be the square of the spatial lattice 
interval. The errors in the time derivatives are therefore of the same order as in the spatial 
derivatives. In order to proceed to much higher resolutions implicit differencing, such as the 
Crank-Nicholson scheme should be used, or the time steps would become prohibitively small. 
However, for the resolutions we have used here, the explicit method is quite manageable. 

B.4 Convergence tests 




Figure 14: (a) /ie(1,— 1) and (b) /ie(0.5,0.5) as functions of the number of lattice points, compared 
to a fit assuming second order scaling to the continumii, fitted using the highest two resolution points. 
Then (c) shows the value of c such that h{0, 0) vanishes, the extrapolated continuum value of which 
is 0.592016. 

In both the complex and symplectic implementations various resolutions were used to 
compute the Kahler-Einstein metric, both to check convergence to the continuum and estimate 
error. 

Taking the example of the symplectic implementation, we uniformly covered the co- 
ordinate square < xi < 1,-1 < X2 < 0. Various lattice sizes were used, including 
25 X 25, 50 X 50, 100 x 100, 200 x 200, 400 x 400. The calculations were performed on a desktop 
computer, with the lowest resolution 25 x 25 taking seconds to run, and the highest resolution 
400 X 400 taking many hours. We estimate /ie as h for sufficient flow time that the update of 
/i in a time step is of order the machine precision. Since second order differencing was used to 
implement the flow locally, we expect that any quantity measured, say O, should scale to the 
continuum as Ocontinuum + A^Ocorrection + ©(A^). Values of the relaxed function h at different 
coordinate locations were used to check this scaling, and indeed give this consistent second 



order scaling behaviour. In Fig. 14 we give an example, plotting the value of the boundary 



point /ie(1,— 1), and also the value of an interior point /ie(0.5, — 0.5) where we note that in 



the flow, h{0,0) is fixed to zero by the appropriate choice of c, which is shown in Fig. 14c. 



Using the two highest resolution points we fit the second order scahng behaviour above, and 
see a very good fit to the lower resolution points, with /ie(1, —1) = —0.2439 + 0.077A^ and 
/ie(0.5,-0.5) = -0.0570 + 0.014A2. These fits indicate the error in the value of h-^ at a 
point is about 10~^ for the highest resolution 400 x 400 grid calculated in the symplectic 
implementation. 
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Figure 15: A plot of the value of the Kahler potential at the center of the hexagon as a function of 
A. The four data points correspond to our 40 x 40, 80 x 80, 160 x 160, and 250 x 250 grids. The linear 
fit was made with the two data points with smallest A, giving 0.296008 + 0.1898A^ 



Likewise in the complex coordinate implementation, various resolutions were computed, 
up to 250 X 250. The value of the potential at the center of the hexagon, = l,r; = 1) is 
plotted in Fig. |l^ and we see that the convergence to the continuum value is quadratic in A, 
again consistent with the second order spatial finite differencing. This and other such tests 
suggest our best 250 x 250 grid in the complex case is accurate to about a part in 10^. 

Note that in both the symplectic and com- x 
plex coordinates, the value of 2/-|-c at the center 
of the hexagon was found to be 0.592016. 

The diffusion flow used to study the eigen- 
functions of the scalar Laplacian was also differ- 
enced to second order accuracy. The flow was o.oi 0.02 0.03 \ 0.04 
simulated using the same explicit method, with 
time step ^A^, and for the same resolutions as 
above. Using the lowest non-constant eigenfunc- 
tion which transforms trivially under the action 
of the U{1)'^ and Dq isometries, we plot in Fig. 16 



-6.32 
-6.322 
-6.324 
-6.326 
-6.328 

-6.332 
-6.334 



the eigenvalue extracted for different resolutions, 
using the same initial data for the diffusion, and 
again fitting second order scaling to the two high- 
est resolutions. We see consistent second order 
scaling behaviour. 



Figure 16: The eigenvalue as a function of 
number of lattice points for the lowest non- 
constant eigenfunction which transforms triv- 
ially under the hexagon and U{l)xU{l) sym- 
metry. Again we see good agreement with a 
second order scaling fit function. 
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C. Dg-invariant polynomials 



In this section, we describe the set of polynomials 

V = J2bn,mX'lx^ (C.l) 

n,m 

invariant under the dihedral group acting on the hexagon. 

Our dihedral group is generated by two elements Ri and R2 (see ( |3.5C1| )). If we take a 
linear combination of the then Ri{v ■ x) = v ■ Ri ■ x = {R\v) ■ x, i.e the Ri act on v 

by their transpose. 

Now, i?2 has eigenvalues e'^*/^ and e""^*/^. A convenient basis of eigenvectors is 

The basis is convenient because a polynomial that is symmetric in ai = vi ■ x and a2 = V2 ■ x 
is symmetric in xi and X2- 
Consider the polynomial 

^ = ^c„,^a>^. (C.2) 

n,m 

Symmetry under interchange of ai and a2 requires Cn,m = Cm,n- Moreover, we have 

i?2Ka^) = exp (^y (n - m)) (C.3) 

which implies that n— m = mod 6. We conclude that the most general polynomial invariant 
under the group action can be decomposed into a sum of polynomials of the form a]^a™+a^a2 
where n — m = mod 6. 

For example, two important invariant polynomials are 

0102 = xf + X\X2 + X2 = U , 

^{af + 4)^ = xjxlixi + X2f = V . 

We now argue that any invariant polynomial can be decomposed into sums and products 
of U and V. Assume a"a2* + a™a2 ^^^^ invariant by R2 and assume m is the minimum of 
m and n. Then 

CL1CL2 + di CL2 — \fl\(l2) j — U (a-|^ j 

One can prove inductively that a""™" + 02""* written in terms of U and V whenever 

(n — m) is a multiple of 6. So U and V generate the polynomials invariant under Ri and R2- 
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D. Convergence and error estimates of the constrained optimization method 



In this appendix we give convergence results for the smooth part of the symplectic potential, 
h{xi,X2) determined in Section |6|. At all points in the hexagon domain the value of h was 
observed to converge quickly with increasing numbers of terms taken in the series expansion 
at the origin. In Fig. ^ we plot the value of /i at 3 points. Note that one of these points 
lies on the boundary of the hexagon. All other points checked, both in the hexagon interior 
and on the boundary, gave qualitatively similar convergence. In the figure we also compare 
the data with the extrapolated continuum results for the Ricci flow given in the previous 
appendix (suitably adjusting for the different value of c). We see excellent agreement. We 
may also estimate that the eighteenth order results and the estimated infinite order result 
differ at the 1 part in 10^ level. 
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Figure 17: (a) h{0,0), (a) /i(0.5, — 0.5) and (b) h{l,—l) as functions of the number terms in the 
expansion. The dashed lines show the extrapolated continuum from the Ricci flow results. 



References 

[1] S.-T. Yau, "Calabi's Conjecture and some new results in algebraic geometry," Proc. Nat. Acad. 
Sci. 74, 1798 (1977); 

S.-T. Yau, "On the Ricci curvature of a compact Kahler manifold and the complex 
Monge-Ampere equations," Comm. Pure App. Math. 31, 339-411 (1978). 

[2] M. Headrick and T. Wiseman, "Numerical Ricci-flat metrics on K3," Class. Quant. Gray. 22, 
4931 (2005) [arXiv:hep-th/0506129]. 

[3] S. Donaldson, "Some numerical results in complex differential geometry," 
arXiv:math.DG/0512625. 

[4] M. R. Douglas, R. L. Karp, S. Lukic and R. Reinbacher, "Numerical Calabi-Yau metrics," 
arXiv:hep-th/0612075. 

[5] M. R. Douglas, R. L. Karp, S. Lukic and R. Reinbacher, "Numerical solution to the hermitian 
Yang-Mills equation on the Fermat quintic," arXiv:hep-th/0606261. 

[6] G. Tian, "On Kahler-Einstein metrics on certain Kahler manifolds with ci{M) > 0," Invent. 
Math. 89 (1987) 225-246; 

G. Tian and S. T. Yau, "On Kahler-Einstein metrics on complex surfaces with Ci > 0," 
Commun. Math. Phys. 112 (1987) 175-203. 



- 40 - 



[7] Y.T. Siu, "The existence of Kahler-Einstein metrics on manifolds with positive anticanonical 
line bundle and a suitable finite symmetry group" Ann. of Math. (2) 127 (1988), no. 3, 585-627. 

[8] B. Feng, A. Hanany and Y. H. He, "Phase structure of D-brane gauge theories and toric 
duality," JHEP 0108, 040 (2001) [arXiv:hep-th/0104259]. 

[9] A. Hanany and A. Iqbal, "Quiver theories from D6-branes via mirror symmetry," JHEP 0204, 
009 (2002) [arXiv:hep-th/0108137]. 

[10] C. E. Bcaslcy and M. R. Plesser, "Toric duality is Seiberg duality," JHEP 0112, 001 (2001) 

[arXiv:liop-th/0109053]. 

[11] I. R. Klebanov and A. A. Tseytlin, "Gravity duals of supersymmetric SU(N) x SU(N+M) gauge 
theories," Nucl. Phys. B 578, 123 (2000) [arXiv:hep-th/0002159]. 

[12] I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: Duality 
cascades and chiSB-resolution of naked singularities," JHEP 0008, 052 (2000) 
[arXiv:hep-th/0007191]. 

[13] M. Abreu, "Kahler geometry of toric manifolds in symplectic coordinates," to appear in Toric 
Varieties in Algebraic Geometry and Physics V. Batyrcv (cd.), AMS, math.DG/0004122 

[14] G. Tian and X. Zhu, "Convergence of Kahler- Ricci Flow," J. Amer. Math. Soc, posted Nov 
17 2006, 8-0894-0347(06)00552-2 (2006) 

[15] K. Hori and A. Kapustin, "Duality of the fermionic 2d black hole and N = 2 Liouville theory as 
mirror symmetry," JHEP 0108, 045 (2001) [arXiv:hep-th/0104202]. 

[16] D. Garfinkle and J. Isenberg, "Critical behavior in Ricci flow," arXiv:math.DG/0306129; 

D. Garfinkle and J. Isenberg, "Numerical studies of the behavior of Ricci flow," in Geometric 
evolution equations, vol. 367 of Contemp. Math., pp. 1-3-114. Amer. Math. Soc, Providence, 
RI, 2005. 

[17] M. Headrick and T. Wiseman, "Ricci flow and black holes," Class. Quant. Grav. 23, 6683 
(2006) [arXiv:hep-th/0606086]. 

[18] M. Headrick and T. Wiseman, to appear 

[19] X. Wang and X. Zhu, "Kahler-Ricci solitons on toric manifolds with positive first Chern Class," 
Adv. Math. 188 (2004) no. 1, 87-103 

[20] V. Baytrev and E. Sclivanova, "Einstcin-Kahler metrics on symmetric toric Fano manifolds," J. 
Reine. Agnew. Math. 512 (1999), 225-236. 

[21] www.stanford.edu/~headrick/dP3/ 

www.phys.washington.edu/'herzog/ dP3 / 

[22] S. Donaldson and J. Keller, to appear 

[23] S. B. Giddings, S. Kachru and J. Polchinski, "Hierarchies from fluxes in string 
compactiflcations," Phys. Rev. D 66, 106006 (2002) [arXiv:hep-th/0105097]. 

[24] S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, "De Sitter vacua in string theory," Phys. 
Rev. D 68, 046005 (2003) [arXiv:hep-th/0301240]. 

[25] S. Kachru, R. Kallosh, A. Linde, J. M. Maldaccna, L. McAllister and S. P. Trivedi, "Towards 
inflation in string theory," JCAP 0310, 013 (2003) [arXiv:hep-th/0308055]. 



- 41 - 



[26] E. J. Copeland, R. C. Myers and J. Polchinski, "Cosmic F- and D-strings," JHEP 0406, 013 
(2004) [arXiv:hep-th/0312067]. 

[27] K. Altmann, "The versal deformation of an isolated toric Gorenstein singularity," 
arXiv:alg-geom/9403004. 

[28] D. Berenstein, C. P. Herzog, P. Ouyang and S. Pinansky, "Supcrsymmetry breaking from a 
Calabi-Yau singularity," JHEP 0509, 084 (2005) [arXiv:liep-tli/0505029]. 

[29] S. Franco, A. Hanany, F. Saad and A. M. Uranga, "Fractional brancs and dynamical 
supcrsymmetry breaking," JHEP 0601, Oil (2006) [arXiv:liep-th/0505040]. 

[30] M. Bertolini, F. Bigazzi and A. L. Cotrone, "Supcrsymmetry breaking at the end of a cascade of 
Seiberg dualities," Phys. Rev. D 72, 061902 (2005) [arXiv:hep-th/0505055]. 

[31] A. Brini and D. Forcella, "Comments on the non-conformal gauge theories dual to Y(p,q) 
manifolds," JHEP 0606, 050 (2006) [arXiv:hep-th/0603245]. 

[32] M. Grana and J. Polchinski, "Supersymmetric three- form flux perturbations on AdS(5)," Phys. 
Rev. D 63, 026001 (2001) [arXiv:hep-th/0009211]. 

[33] Y. Matsushima, "Sur la structure du groupe d'homeomorphismes analytiques d'une certaine 
varietc kaehlerienne," Nagoya Math. J. 11 (1957), 145-150. 

[34] S. Franco, Y. H. He, C. Herzog and J. Walcher, "Chaotic duality in string theory," Phys. Rev. 
D 70, 046006 (2004) [arXiv:hep-th/0402120]; 

S. Franco, Y. H. He, C. Herzog and J. Walcher, "Chaotic cascades for D-branes on 
singularities," arXiv:hep-th/0412207. 

[35] C. P. Herzog, Q. J. Ejaz and 1. R. Klebanov, "Cascading RG flows from new Sasaki-Einstein 
manifolds," JHEP 0502, 009 (2005) [arXiv:hep-th/0412193]. 

[36] C. P. Herzog and I. R. Klebanov, "On string tensions in supersymmetric SU(M) gauge theory," 
Phys. Lett. B 526, 388 (2002) [arXiv:hep-th/0111078]. 

[37] V. Guillemin, "Kahler structures on toric varieties," J. Diff. Geom. 40 (1994) 285-309. 

[38] W. Fulton, Introduction to Toric Varieties, Princeton University Press, 1993, p 92. 



- 42 - 



